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Mach number 

pressure 
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radius of curvature 
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temperature 
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density 

angle of the flow to free -stream direction 

Subscripts 

immediately behind shock 

cone surface 
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effective value 

swept -cylinder stagnation line 

free stream 

normal component to shock 
maximum value 
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APPROXIMATION FOR DISTRIBUTION OF FLOW 


PROPERTIES IN THE ANGLE -OF -ATTACK 
PLANE OF CONICAL FLOWS 

Joseph W. Cleary 
Ames Research Center 


SUMMARY 


An approximation is made for the hodographs of stagnation- line 
streamlines in the angle-of -attach plane of conical flows that yields explicit 
equations for distributions of velocity and flow angle. It is assumed that 
shock angle and surface inclination are known and the flow is inviscid and 
isentropic. In essence , the approximation replaces streamline hodographs 
of conical flows by circular arcs. Streamline hodographs of wedge and swept- 
cylinder flows are obtained from the solution as special limiting cases of 
conical flow. The approximation introduces no limitation on angle of attack 
or body shape as long as shocks are straight and flows are conical. 

Comparisons of the approximation with exact numerical solutions for the 
flow over right -circular cones demonstrate agreement for Mach numbers from 
one to infinity over a wide range of cone angles and angles of attack. 
Comparisons with experimental results show applicability of the approximation 
to conical flows other than those for right -circular cones. 


INTRODUCTION 


Theoretical aspects of conical flows are of fundamental interest to 
aerodynamic research and to the more practical problems of vehicle design. 

The direct application to lifting entry configurations has enlivened interest 
in conical-flow theory for a wide range of angles of attack and body shapes. 
The flow In the angle -of -attack plane is of unique interest in conical flow 
applications since flap-type controls may be placed on the windward stagna- 
tion line in order to achieve greater effectiveness from compression of the 
flow. Moreover, heating is apt to be most severe on the windward stagnation 
line . 


Knowledge of the compression process in the angle -of -attack plane is 
required to adequately estimate control effectiveness and heating. Although 
inviscid axisymmetric solutions by numerical methods are available (e.g., 
refs. 1 to 3) for estimating the compression process at 0° angle of attack, 
only limited solutions have been published for estimating the compression of 
the flow at angle of attack (refs, k and 5). Moreover, at present, numerical 
solutions of the flow generally require that three-dimensional bodies be 
smooth without sharp edges. Since conical-type bodies with sharp edges on 
the windward surface may prove useful for achieving high drag at high angles 



of attack during entry, some method is needed for calculating the compression 
process of conical flows in general. 

The purpose of the present study is to present an approximate solution 
for the compression process of the flow in the angle -of -attack plane. A 
semiempirical approach is used to develop closed-form equations for the 
distribution of flow properties between the shock and the body in this plane. 
While it is assumed that surface inclination and shock angle are known, no 
limitation is made on angle of attack or body shape to the extent that shock- 
wave elements are straight and the flow is conical. Comparisons are made 
with numerical solutions and with experiment to demonstrate accuracy and 
general applicability. 


ANALYSIS 


Description of the Approximation 
and Presentation of Equations 

The intent of the present analysis is to develop an approximation for 
the distribution of properties of the flow in the angle -of -attack plane of 
conical flows. It is assumed at the outset that the shock angle 0 2 and 
surface inclination u> c are known and that, as depicted in figure 1, the 
flow between the shock and the body in the angle -of -attack plane is inviscid 
and isentropic. Before the approximation is described, the true relationship 
between velocity vectors and angles pertinent to conical flows will be exam- 
ined. This relationship can be shown most clearly by a plot of the velocity 
component v as a function of u (hodograph coordinates). Such a plot can 
also be considered a polar diagram with radius vector V given as a function 
of the turning angle of the flow u>. It is expedient to consider the hodo- 
graphs of axisymmetric flows first and then consideration will be given to 
other conical flows. 

Axisymmetric conical flows - Figure 2 shows the hodograph of a streamline 
of a typical axisymmetric conical flow. In figure 2 , R is the local radius 
of curvature of the streamline hodograph at an arbitrary point p in the 
conical field. Note that since R is normal to the streamline hodograph, 
an extension of R makes the conical coordinate angle Q with the free- 
stream direction. Furthermore, points 2 and c correspond to conditions at 
the shock and cone surface, respectively; points 2 f , p 1 , and c f represent 
end points of R that correspond with points 2 , p, and c, respectively. 

In the present case the magnitude of R was evaluated by an exact equation 
given by Busemann (refs. 6 and 7 ). Busemann*s equation was used to determine 
R and the locus of end points of R from an exact numerical solution of 
the flow. While the differential equation of the streamline hodograph is 
known (ref. 7 ), no general solutions are known. Therefore, an approximation 
is considered appropriate. 

In essence, the present approximation replaces the exact hodograph by a 
circular arc. Details of the procedure are shown in figure 3 where R is 
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replaced by a constant effective radius of curvature R whose center is at 
H and K. The circular arc passes through point 2 and has its center at the 
intersection of R 2 and R c (or an extension of V c ). Since R and R are 
both normal to the streamline hodograph at point 2, the exact velocity and 
slope dv/du are retained in the approximation at the shock wave. While the 
exact value of dv/du at the cone surface is also retained, the surface 
velocity is in error by an unknown amount. 

The circular-arc approximation is given by equation (l). 

(u - H) 2 + (v - K) 2 = R 2 (l) 


Equations for R, H, and K are derived in appendix A and are given by ( 2 ), 
(3), and ( k ), re spe ct i ve ly : 


sin(u c - uj 2 )cos 0 2 

R = v ro (2) 

sin(0 2 - co c )cos(0 2 - to 2 ) 


tan(0 2 ~ ^ 2 ) 

H = V * cos 0 2 cos to (3) 

OO / \ V 

sin(0 2 - to c ) 


tan(0 2 " w s) 

K = — cos 0 2 sin to c (4) 

sin(0 2 " w c) 


where 0 2 , <o 2 , and to c = 6 C are known quantities. Furthermore, it is shown 
in appendix A that b^; substitution of u - V cos <0 and v = V sin w and the 
known values of R, H, and K in equation (l), V/Voo is given as a function 
of (to c - to) by equation ( 5 ). 



cos 0 2 tan(0 2 " t^) 
sin(0 2 - ^c) 


(cos(co c 


— 


, \ 



1/2 

sin 2 (to c 

- U) 2 ) 


COS 2 ((O c - to) - 1 + 



sin 2 (0 2 

- w 2 ) 






(5) 


Since additional pertinent equations are derived in appendix A, it is suffi- 
cient here to summarize the final results. 
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Velocity distributions normalized by Voo and Vc are given as functions 
of the conical coordinate ( 0 - w c ) by equations (6) and (1) , respectively. 


— = cos 0 2 tan 0 2 tan ~ (d 2 -u c ) + 1 jl + 


2B 


(l-B)‘ 


[i - cos(e - u )] 


1/2 


tf- = 1 + 


2B 


(1 - b) 


[i - cos(e - u c )] 


1/2 


(6) 

( 7 ) 


where the parameter B is given by equation (8) 

/ \ sin(e 2 - to c ) 

b = oos(e 2 - » c ) - (p j P2 ) tan e - 


( 8 ) 


and the density ratio across the shock P^/Pg by equation (9) (see ref- 8) 



( 9 ) 


Properties of the flow at the cone surface are of particular interest and can 
be evaluated from equation (6) by specifying that 0 = co c to give 
equation (10) 


tt- = cos 0 2 
*00 


P 1 
— tan 0 2 tan — (0 2 - w ) + 1 


( 10 ) 


The local flow angle (to c - w) is given by equation (ll) as a function of the 
conical coordinate ( 0 - w c ). 


io c - to = tan 


■ l B sin(0 - mo) 

1 - B cos(0 - w c ) 


( 11 ) 


An interpretation of how to apply equations (l) to (ll) to conical flows in 
general will now be given. 

Other conical flows - Since in the angle- of-at tack plane conical flow is 
characteristically isentropic, there is a close relationship between the 
axisymmetric solution and conical flows in general. Therefore, it is assumed 
in like manner that a circular arc approximates the streamline hodograph of 
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the general conical flow. Within this assumption equations (l) to (ll) are 
applied to the general case a > 0° by assigning w c = 5 C -i- a and 
w c = & G - a on the windward and leeward stagnation lines, respectively. If 
on the leeward side the singular line, for which entropy is multivalued, has 
lifted from the surface, then oj c = v - a where v is the angle between the 
cone reference axis and the singular line. 

The usefulness of the circular-arc assumption for a > 0° will be 
scrutinized later when the approximation is compared with exact numerical 
solutions. At present, support for it can be given by examining two limiting 
conical flows in the hodograph plane; these are the flow over a wedge and the 
flow in the sweep plane of swept cylinders. Figure k shows the holographs of 
these flows and of their streamlines together with that of a typical conical 
flow. The hodograph for wedge flows is, of course, the shock polar while 
that for swept cylinder flows is the trace of the component of velocity 
parallel to the wave and is given by the semicircle Vg/v^ = cos 0 2 . 

For a specified 0 2 , the streamline hodograph for wedge flow is point 2 
(u and v constant) on the shock polar; that for the swept cylinder is a 
straight-line extension of the change in velocity through the shock wave from 
point 2 to point s on the swept- cylinder polar (see fig- 4 ). Results from 
a limited number of exact numerical solutions for conical flows (refs. 1 to 5) 
demonstrate that V c must terminate within the region bounded by the shock 
and swept -cylinder polar s and therefore, this is accepted as a factual con- 
dition. Moreover, numerical solutions show that for a specified 0 2 and oo, 
the velocity of conical flow exceeds that for a swept cylinder, but this 
velocity is, of course, less than that for a wedge having the same 0 2 . 
Therefore, it can be shown that the inequality given by (l2) applies (see 
fig* ^)* 


V s < V 


cos 9^ 


00 cos (0 2 - w c ) 


< < V < V. 


( 12 ) 


For constant 0 2 it is apparent that there Is a spectrum of conical 
flows whose limits are the flow over a wedge and a swept cylinder. For these 
flows equation (10) predicts a single distribution of ^c/^oo a function 
of <a> c . Wedge and swept -cylinder flows are limiting cases for which <a) c = oo 2 
and w c = e 2 , respectively, as can be verified by substitution in equation (10) 
to yield the respective exact results V c /v = - cos 0 2 /cos(0 2 - w 2 ) 

and V c /Voo = Vg/Yoo - cos @ 2 . Similarly, as demonstrated in appendix B, it can 
be shown that distributions of v/v^ given by equation ( 5 ) contain wedge and 
swept -cylinder distributions as special cases. 

Finally, consider the effects of and 7 on the present approx ima tion . 
As Mx> increases or 7 decreases, the shock polar grows in size and in the 
limit = 00 and 7 = 1 coincides with the swept -cylinder polar. It is 
evident that for this limit, w c = and V c = V s = V 2 = cos 0 2 . 
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Furthermore, as the limit is approached the length of the streamline hodograph 
of V is shortened and V c V 2 . Consequently, the accuracy of the approxi- 
mation is improved, as will he demonstrated later hy comparisons with exact 
numerical solutions. First, however, before demonstrating accuracy, equations 
will be given for other flow properties in a form to expedite estimates. 

Equations of stat e prop erties- Since the flow in the angle- of -attack 
plane is isentropic, other flow properties can be estimated locally by use of 
equation ( 6 ) in conjunction with Bernoulli’s equation and equations from 
supersonic flow and oblique-shock theories. Equations for other flow proper- 
ties were derived from equations given in reference 8 and are summarized below 
for completeness. 


where 


1 ± 1 is _ q 

T_ 7 - 1 P g 2 

Too y t i Ppp 

7 - 1 P 2 



1 

p > 2 


7+1 P oo 

7 - 1 Pa 

7+1 Pcs _ 
7 - 1 P 2 



l 

7-1 


_p_ jt 

P p T 

-^00 CO £> 



(13) 

(14) 

(15) 

(16) 


To estimate Q note that [(v/V 00 )/cos & 2 ] can be obtained directly from equa- 
tion ( 6 ) or at the cone surface from equation (lO) . Since it can be shown 
that at the surface = ( Q^J P 2) 2 and_ 0 _ for wedge and swept -cylinder flows, 

respectively, it is evident that 0 < Q < 1 for all conical flows. 

Since real-gas flows are of practical interest, the application of 
equations ( 13 ) to (l 6 ) to the specific case of a real gas in thermodynamic 
equilibrium will be considered briefly. For simplicity, it is assumed that 
7 is constant, since in an actual case the variation of 7 along a stream- 
line, behind the shock, is small. An effective value of 7 can be estimated 
from known free- stream conditions and the density ratio across the shock and 
is given by equation (l 7 )« 
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( 17 ) 


( 2 A€n ) - M-g * j 

(p./p*) - 1 


With known, flow properties are estimated from equations (l 3 ) to (l6) 

by replacing 7 with its effective value from (17). 

RESULTS AND DISCUSSION 


To demonstrate accuracy of the present approximation, the predicted 
results will be compared with exact numerical solutions for conical flows. 
Surface properties will be considered first; then the predicted distributions 
of V and oo will be compared. Since all surface properties of the flow may 
be of interest, the relative errors between them for a specified error in 
velocity may be pertinent and will be considered briefly. Starting with the 
isentropic flow relation between p c /p 2 and V c given by equation (18) 


7 

2 - V 2 \ 7-1 

'Tv) 

it can be shown that the relative error between P c /p 2 and V c may be estimated 
by means of equation (19). 




d ( p c / p g)/(p c / p g) 

dV /V 
c c 


(27/7 - i)(v c /v m ) g 
1 - (v /v ) 2 

x c nr 


(19) 


For a specified error in V c , equation ( 19 ) predic ts that the e rror in p c /p 2 
■will num erically excee d that in V c if V c /v m > J (j-l) / ( 37 - l) whereas if 
V c /V m < J (7-1)/ ( 37 _ l) the converse is true. Moreover, it can be shown that 
errors in p c /p 2 are greater than in p c /p 2 , T c /T 2 , and c c /c 2 in the respec- 
tive proportions ( 2 7/7 - l) to (2/7 - l), 2 , and 1. Therefore, a comparison 
of predicted V c and p c with exact solutions should suffice to give an upper 
bound on errors in most flow properties of interest. This comparison will be 
made on a limited basis since only solutions for conical flows generated by 
right -circular cones are available in a form suitably accurate for comparison. 
Later other conical flows will be compared with experiment. 
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Comparison of Predicted Surface Properties 
With Inviscid Numerical Solutions 

Results from the approximation will he compared initially with 
axi symmetric solutions and then comparisons will he made for a > 0°. 

Comparisons with axisymmetric s olutions - Surface velocity and pressure 
will he compared with results from solutions given hy references 1, 2, and 9. 
References 1 and 2 give tabulated results from numerical solutions of the 
Taylor-Maccoll equation for conical flow, while reference 9 provides solutions 
hy the method of characteristics. The present approximation is compared with 
solutions hy Sims (ref. l) for 7 = 1.4 in figure 5. Here the approximations 
for V c and p c are ratioed to exact values, and constant values of these 
ratios are shown as contours of accuracy for 0 < S c < 30° and 1.5 < M < 20. 

In like manner comparisons are shown in figure 6 with solutions hy Kopal 
(ref. 2) for 7 = 1.405 and 1.333* JIhese comparisons are made for 
5° < & c < 50° hut are limited to 1 < < 12 because the number of solutions 

given was insufficient to define accuracy contours at greater M^. 

Figures 5 and 6 demonstrate that accuracy of the approximation Improves 
significantly with increasing M^. In general for M w > 5 and b Q < 5O 0 veloc- 
ity is predicted within 0.1 percent and, except for & c 5° * pressure is 
predicted within 1 percent. While accuracy is minimum for small cone angles 
(5 C < 5°) inaccuracies accrue also as 6 C approaches shock wave detachment. 
Rigorous comparisons at detachment could not he made since available numerical 
solutions are given graphically; however, comparison with the method of refer- 
ence 10 for 3y^ = 00 and I.67 < 7 < 1.2 indicate agreement with present predic- 
tions of velocity and pressure to about four significant figures. The present 
approximation appears least accurate as 1. From a comparison with 

results from reference 1 for Mq = 1 figure 7 shows that predicted velocity 
and pressure have minimum accuracy at = 1.06 where both are in error about 
5 percent. 

Figure 6 shows that decreasing the specific-heat ratio improves accuracy 
significantly except when is small. The effects of specific-heat ratio 

were investigated in more detail and are shown in figure 8 where comparisons 
are made with solutions hy the method of characteristics (ref. 9)* Figure 8 
shows that accuracy of the approximation improves as 7 1 for a wide range 

of cone angles. The applicability of the approximation to flows of a real gas 
in equilibrium for a range of 7^ is demonstrated in figure 9. Here the 
pressure ratio p c /p 2 is compared with -exact solutions hy the method of 
reference 9- For the specified free-stream conditions, figure 9 indicates 
that predicted values of p c /p 2 agree with exact solutions to within 
0.5 percent. 

Comparisons with solutions for a > 0°- A critical assessment of how 
well the approximation applies to the more general case a > 0° can he made 
hy a comparison of predicted distributions of the velocity component e c 
with exact numerical solutions of conical flow from reference 4. As shown in 
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figure 4 e c is the incremental value that V c exceeds the local velocity 
for swept -cylinder flow having the same 9 Z and local to = to c and is given 
by equation (20). 


e c = v c 


v cos ..Qa 

00 cos(0 2 - to„) 


] 


( 20 ) 


Results from cross plotting numerical solutions with constant but 

•varying S c and a. are compared with, the approximation in figure 10. Here 
€ c /y oo is presented as a function of tan(02 - co c ) /tan(0£ “ ^ 2 ) as an abscissa. 
From figure 4 it can be seen that this abscissa is the fraction of the veloc- 
ity component from point s to 2 to where V c crosses this component. 

Figure 10 shows that exact distributions of agree well with 

predicted distributions and moreover for specified tend to give a single 

distribution as predicted by the approximation. The condition for s c /V 00 to 
be maximum is predictable from the present approximation by equating 
de c /d[tan (02 - w c )/ta n(0^ - u^)] to zero and is given by equation (21). 


tan* 


■<* - [ £& - ^>] 4 + 2 j 


tan(0s 




5 ) 


tan “ w s) 


1 = 0 


( 21 ) 


Figure 10 shows good agreement between the predicted maximum and exact 
solutions. 

The approximations for V c and p c are compared with conical-flow 
solutions for a > 0° from reference 4 in figure 11. Velocity and pressure 
are normalized by exact values and are shown as functions of the inclination 
of the stagnation line oo c . Comparisons are made for both windward and lee- 
ward stagnation lines, and since for the windward side o) c = S c t a while 
for the leeward side o) c = 6 C - a, the demarcation between these two cases is 
apparent in figure 11 when a) c = 6 C . Trends of accuracy indicated in figure 11 
appear consistent with those previously demonstrated for axisymmetric flow, 
and generally good agreement of predicted and exact results is demonstrated. 

A more practical demonstration of the usefulness of the approximation is 
afforded by the direct comparison of predicted pressure ratio Pc/p 2 with 
exact values shown in figure 12. Also shown are pressure ratios for wedge 
and swept -cylinder flows which are obtained from the approximation as special 
limiting cases. It is evident from figure 12 that for a > 0° the 
approximation provides good accuracy for a wide range of conical flows. 
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Comparisons of Predicted Distributions of Velocity 
and Plow Angle With Inviscid Solutions 

Axisymmetric distributions of velocity and flow angle given by 
equations ( 7 ) and ( ll) , respectively, are compared with exact solutions in 
figure 13. Likewise for a > 0°, distributions in the angle -of -attack plane 
are compared in figure 14. From figure 13 it can be seen that, in accordance 
with the predictions of surface properties, the predicted distributions of 
velocity and flow angle are least accurate at small S c . In general, if sur- 
face properties are predicted accurately then distributions will also be 
accurate. It should be noted that while figure 13 shows significant differ- 
ences between predicted and exact distributions for small 6 C the actual 
differences in flow properties are not large. For example, the greatest 
difference in flow angle between predicted and exact values shown by the 
distributions of figure 13 is only about 1° for 5 C = 2.5°« 


Comparisons of the Approximation With Experiment 

Comparisons with experiment are, of necessity, limited to investigations 
from which both ©2 and p c (or some other flow property) can be evaluated. 
Accordingly, results from references 11 to 14 have been selected for compari- 
sons with the approximation. These investigations encompass a range of Mach 
numbers from about 5 to 15* It is clear from previous inviscid comparisons 
that for circular cones at small a at least, any significant differences 
between the approximation and experiment for this range of M^ must be 
attributed mainly to viscous effects and inaccuracies in the evaluation of 
0 2 and p c rather than to the small errors inherent in the approximation. 
Except for results from reference 11, shock angles were measured from shadow- 
graphs of the flow; shock angles from reference 11 were estimated from pitot 
traverses of shock layers. These traverses also afford an estimate of 
boundary-layer thickness; therefore, the effects of boundary-layer displace- 
ment were accounted for approximately by measuring surface inclination oj c 
to the boundary-layer edge, and to the shear line of the lee separated flow. 

On the other hand, results from references 12 to 14 are compared on the basis 
of w c measured from the cone surface and are denoted inviscid. In addition, 
for these cases w c was arbitrarily increased 1° to show the effect of 
boundary layer displacement. For all cases, pressure was measured at the cone 
surface and the pressure ratio P c /P 2 has, therefore, been selected as a 
basis for comparison. 

In figure 15 Pc/P 2 predicted by the approximation is compared with 
results from references 11 and 12 for a 15° semiapex cone. The agreement 
between experiment and the approximation is reasonable for both air and helium 
on the windward surface; but on the leeward side, greater differences are indi- 
cated for air (fig. 15 (a)), probably because of difficulty in accurately mea- 
suring p c . A similar comparison is shown in figure 1 6 with experimental 
results from reference 13 for S c = 5° • While the approximation appears to 
generally underestimate Pc/P 2 > reasonably good agreement is indicated for 
5° < w c < 45°. 
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Comparisons for an elliptic cone for two different ratios of major -to - 
minor axes are shown in figure 1J from results given in reference 14. 
Although the approximation underestimates experimental values slightly, good 
agreement is indicated for a wide range of oo c . The approximation appears 
reasonably valid even at large w c where 0^ has the strong shock solution 
and the flow in the shock layer is subsonic and therefore affected by the 
model base. Since the shock waves appear straight in the shadowgraphs of 
reference 14, even when Me < 1, it is apparent that the flow continues to be 
nearly conical even though model length is finite. Moreover, shadowgraphs of 
a 15 ° semiapex cone in helium at ^ =* 15 shown in figure 18 demonstrate that 
the shocks tend to remain essentially straight even though 0 2 90 °. It 

appears from these results that if the body is sufficiently slender the flow 
remains nearly conical near the apex even at large a and the approximation 
can be applied. 


CONCLUDING REMARKS 


An approximation is developed for the distribution of flow properties in 
the angle -of -attack plane of conical flows. It is assumed that shock angle 
and surface inclination are known and the flow is invlscid and isentropic. 
Comparisons of the approximation with exact solutions and experiment demon- 
strate applicability for a wide range of conical flows. Equations are 
presented from which estimates of flow properties can be made. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, April 23 , 1970 
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APPENDIX A 


DERIVATION OF EQUATIONS FOR DISTRIBUTIONS 
OF VELOCITY AND FLOW ANGLE 

PROPERTIES OF THE CIRCULAR ARC APPROXIMATION 


The purpose of this appendix is to derive explicit equations for 
distributions of velocity and flow angle by approximating the streamline hodo- 
graphs of axisymmetric conical flows by circular arcs. Pertinent velocity vec- 
tors and angles are shown in sketch (a) . The center^of the circular arc_with 

radius R is at the point (H, K) 
which is the intersection of an 
extension of V c and a normal to 
the component 2-s at point 2. By 
constructing a normal to V c that 
passes through point 2 and by 
resolving components it is evident 
that 

R sin (©2 - w c) = ^2 sin(oo c - oj^) 

(Al) 

and since 



Sketch (a) 


V 2 = V w 


cos 9 g 

cos(02 " w 2 ) 


(A2) 


it follows that 


r - v M 


cos 9g sin(co c - a) g ) 

cos( 0 2 - W2)sin(02 - w c) 


(A3) 


In addition 


H = V 2 cos 


+ R co s 0 2 


(Alt) 


K = V 2 sin oj 2 + H sin 0 2 


(A5) 


Substituting (A2) and (A3) in both (Ak) and (A5) and employing known 
trigonometric identities yields (A 6) and (A 7 ), respectively. 


H 


tan(02 - 0 J 2 ) 
sin (02 " w c) 


cos COS 


co c 


(A6) 
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cos 0 2 sin u> c 


(AT) 


Velocity as a Function of (u) c - w) 

The equation of the circular arc is given hy (AS). 

(u - H) 2 + (v - K) 2 = R 2 (A8) 

Substitution of u = V cos a) and v = V sin w in equation (A8) yields the 
quadratic equation in V given by (A9) 

V 2 - 2(5 cos a) + K sin w)V + H 2 + K 2 - R 2 = 0 (A9) 

Substituting equations (A3), (A6), and (A7) for R, H, and K, respectively, in 
equation (A9), solving the quadratic, and employing trigonometric identities 
yields V/Voo as a function of (o) c - w) and this is given by equation (A 10) . 



(A10) 

It is apparent that the minus sign before the radical is the appropriate 
choice here since the plus sign yields velocity to the opposite side of the 
circle . 

Flow Angle as a Function of (0 - o) c ) 


By erecting normals through the point H and K to extensions of V and 
and resolving components it can be shown that equation (All) applies 


sin(0 2 - u) 2 )sin(io c - to) = sin(u) c - ufe)sin(0 - w) 


Equation (All) can be recast 1 as (A12) where (u c - oo) is an explicit function 
of the conical coordinate (0 - oo c ) 

f \ 1 - B COS(0 - U) G ) { k \ 

COS(03 C - U>) = ^ TJZ ( A12 ) 

[B 2 + 1 - 2B cos(0 - o> c )] 


■^Rote that sin(0 - oo) - sin[ (0 - w c ) + (o) c - co)], 
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and where 


sin(a)c - (*>g) 
sin(0 2 “ w 2 ) 


A simpler alternate for (A12) is given by (Al4) 




0) - 


tan 


B sin(0 ■ a y») 

1 - B cos(0 - co c ) 


Velocity as a Function of (0 - co c ) 


(A 13) 


(AU) 


Substitution of (A 12) in (A10) yields V/Voo as a function of the conical 
coordinate (0 - u) c ) and this is given by equation (A15) 


cos 0 2 Jtan(0 2 ~ w 2 )tan (0 2 - w c )+ lj|l + -p— [1 - cos(0 - co c )]J 1/2 

(A15) 

Using oblique -shock theory in conjunction with an alternate equation for B, 
equation (A15) can be recast independent of to 2 . From oblique-shock equations 
given in reference 8 it can be shown that equation (Al6) applies 2 


r-' 

tan (0 2 - w 2 ) = q 22 tan 0 2 


(A16) 


where 




7 + 1 Un + ^ 


(A17) 


A combination of equation (Al6) with an expansion of (A13) yields equa- 
tion (A 18) which is independent of u) 2 . 


B = cos(0 2 - o> c ) 


sin(0 2 - w c ) 

rpoo/p 2 ) ta n 0 2 


(A18) 


Substitution of (Al6) in (A15) with B given by equation (Al8) yields equa- 
tion (A 19) which is independent of u^. 3 


2 From equation (Al6) and figure 4 it is evident that p^/p^ equals the 
ratio of the velocity component s - 2 to the free -stream component normal to 
the shock wave. 

^Equation (Al4) is made independent of w 2 by using (Al8) for B. 
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V o 

— = cos 0 2 

V oo 


p!f tan Q s tan | (S 2 - w c ) + 1 ( 1 + JJ ? B ) g [1 - cos(0 - w c )]j 1/2 

(A19) 


At the cone surface 9 - u) Q and equation (A19) reduces to (A20). 

v c r Poo L 1 

— = cos 0 £ ~ tan 0 2 tan - (0 2 - a> c ) + 1 


(A20) 


Equation (A19) when divided by equation (A20) yields velocity distribution 
normalized by V c and this is given by equation (A2l) 


X 

V c 



2B 

(1 - B ) 2 


[1 


cos (0 


^c) 



(A21) 
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APPENDIX B 


COMPATIBILITY OF THE APPROXIMATION WITH 
WEDGE AND SWEPT -CYLINDER FLOWS 


Equation (Bl) (identically eg. ( 5 )) will be examined to verify that 
distributions of V/V^ for wedge and swept -cylinder flows are contained as 
special cases. 


V cos 0 2 tan(0p - Wp) 

Voo " sin(0 2 - w c ) 


cos( Wc - w) - [cos £ (w c - w) - 1 + J )'] } 

(Bl) 


For wedge flow w„ = w 2 = u and substitution in equation (Bl) gives 
V/Voo = cos 0 2 /cos(0 2 - ufe) = V 2 /V 00 which is the exact result (see fig. 4). 


For swept cylinders w c = 0 2 and substitution in equation (Bl) gives an 
indeterminate form since the outer bracketed quantity and the denominator are 
both zero. An evaluation can be made by allowing 


F = sin(0 2 - Wc) 


(B2) 


and 


f = cos(w c - w) - [cos 2 (w c - w) - 1 + f— 0 Jgj ] ■ ^ B3) 


and evaluating 


1 . df /dw c 

^ dF/dWc 
w c -^0 2 ' G 

Performing the indicated operation gives 


(b4) 


lim 

w c -»0 2 


_ii f -cos(e 2 - hte) 


cos( 0 2 ~ w c) Lcos(0 2 “ w)sin(02 “ ^)J Lcos(02 “ w)sin($2 “ ^2) J 


H 


c . ° . s ( e £ - ti a) 


(B5) 


and 


~ = cos 0 2 tan(0 2 - wg) I" c° s .( 0 2 ~ ~ 7 ^ — r"| (B6) 

Voo L co s (0 2 - w)sm(0 2 - w 2 ) J 
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Equation (B 6) can Be reduced to 


_V cos 9? 

Voo “ cos(0 2 - w) 

which can be recognized as the exact expression for swept -cylinder flow from 
figure 4. 
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Figure 1.- Example of imriscid conical flow. 



Figure 2.- Streamline hodograph of a typical axisymmetric conical flow 
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n for the streamline hodograph. 



of typical conical flows 






(b) Pressure ratio. 

Figure 5.- Accuracy contours of the approximation from comparisons with exact 

solutions of reference 1; a = 0°. 



(a) Velocity ratio. 


8c. de 9 



0 2 4 6 8 10 12 


Mco 

(To ) Pressure ratio. 

Figure 6.- Accuracy contours of the approximation from comparisons with exact 

solutions of reference 2; a = 0°. 
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(a) Velocity ratio. 



(D) Pressure ratio. 

Figure 7.- Comparison of the approximation with exact solutions from 

reference 1 with ML = 1. nr. = 0°. 








8 C , deg 


Figure 9«“ Comparison of the approximation with real gas solutions for air in 
thermodynamic equilibrium; M = 20 ? V = 23,414- ft/sec, P = 1 lb/sq. ft. 




Approximation' 


Eq (20) 

Eq (21) 


Exact solutions, flagged symbols 
denote extrapolations. 

• Axisymmetric 
a> 0 deg 


S c , deg 

A 10 O 30 

O 15 A 35 

□20 □ 40 

V 25 O 45 




0 .2 .4 .6 .8 1.0 


Tan (0 2 -oj c ) /Tan (0 2 -w 2 ) 


(a) e s = 22. 5 C 


(b) 0 2 = 48° 


Figure 10.- Comparison of approximate incremental velocity with exact 

solutions from reference 4. 
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(a) Velocity ratio . 






(t>) Pressure ratio. 

Figure 11.- Comparison of the approximation with exact solutions from reference 4; a > O' 











P c /P2 



(a) M ra = 10 . 6 ; 7 = 1.4-, reference 11 . 



(b) = lh.9; 7 = 1*67, reference 12. 

Figure 15*“ Comparison of approximate pressure ratio with experimental values 

for a 15° semiapex cone. 






a) c = (B c + a), deg 


Figure 1 6 .- Comparison of approximate pressure ratio with experimental values 
for a 5° semiapex cone; M = 5*03, 7 = 1»4, reference 13 • 




(a) a/b = 0.667; S c = 15 . 2 ° 

P c /P 2 



(b) a/b = 1 . 5 ; S c = 10.3° 


Figure 17.- Comparison of approximate pressure ratio with experimental -values 
for elliptical cones; M = 10, 7 = 1.4, reference 14. 
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